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ABSTRACT 



We report the detection of carbon monoxide (CO) emission from the young 
supernova remnant Cassiopeia A (Cas A) at wavelengths corresponding to the 
fundamental vibrational mode at 4.65 fim. We obtained AKARI Iniraxed Camera 
spectra towards 4 positions which unambiguously reveal the broad characteristic 
CO ro-vibrational band profile. The observed positions include unshocked ejecta 
at the center, indicating that CO molecules form in the ejecta at an early phase. 
'p_j We extracted a dozen spectra across Cas A along the long 1' slits, and compared 

these to simple CO emission models in Local Thermodynamic Equilibrium to 
obtain first-order estimates of the excitation temperatures and CO masses in- 
volved. Our observations suggest that significant amounts of carbon may have 
been locked up in CO since the explosion 330 years ago. Surprisingly, CO has not 
been efficiently destroyed by reactions with ionized He or the energetic electrons 
created by the decay of the radiative nuclei. Our CO detection thus implies 
that less carbon is available to form carbonaceous dust in supernovae than is 
currently thought and that molecular gas could lock up a significant amount of 
heavy elements in supernova ejecta. 

Subject headings: molecular processes - ISM:molecules - supernova remnants:Cas 
A - ISM:dust 
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Introduction 



Molecular emission from supernovae (SNe) remains rarely detected, as the formation 
of even the most stable molecules is difficult in the SN ejecta. The high UV radiation due 
to recombinations, and the energetic electrons from the radioactive decay of ^^Ni and ^^Co 
dissociate molecules efficiently. The presence of molecules in ejecta could thus imply an 
inhomogeneous distribution of material within the ejecta and c ould constrain the degree of 
mixing. Since CO is an effective coolant (ILiu fc Dalgarndll995l ). it enhances dust formation 
and thus the formation efficiency of CO critically affects the dust formation efficiency. Be- 
cause of the rapid decrease in density and temperature during condensation in the ejecta, 
dust formation happens through non-equilibrium condensation. There are then two impor- 
tant key questions for dust condensation in supernovae. First, to determine which molecular 
species are relevant in non-equilibrium condensation and second, whether the gas in the 
ejecta is able to mix with the molecular species. 

There are only about a dozen SNe known where CO has been detected. SN1987A is the 
only SN with spectro s copic detections of both the CO fundamental and first overtone bands 
( IWooden et al.l Il993t ISpyromilio et al.l 119881 ). Other detections of the fundamental band 
have been made with the Spitzer Space Telescope, but these observations cover only the tail 
of CO feature at w avelengths longwards of 5.3 /im for SN2004dj, SN 2005aj, and SN2006et 
(IKotak et al.ll2006l ). More detections have been reported in the first overtone band in the 
near-infrared, wi th observations of SN 1995ad, SN 1998S, SN 1998dl, SN 1999em, SN 1999gi, 
and SN 2000ew (IGerardy et al.ll2002l ). In all cases, the CO detections correspond to times 
between 97- 200 days after the explosion. 

We previously detected CO i n Cas A in near- infrared narrow-band images centered at the 
2.29 /im first overtone CO band (IRho et al.ll2009l ). In this paper, we report the spectroscopic 
detection of the fundamental CO band in Cas A with the AKARI Infrared Camera. 



Observations 



We performed AKARI Infrared Camera (IRC; Onaka et al. 2007) observations towards 
4 positions in Cas A. Three positions are the brightest knots in the near-infrared 2.29 fim 
narrow-band Palomar image (IRho et al.ll2009l ): two positions are towards the northern shell 
(filenames are "nl" and "n2"), one towards the south ("si"), and an addition al position 
towards the central position where unshocked silicon ejecta (sipeak) are detected (IRho et al. 

mi- 



The observations took place in January and July 2009 during AKARI phase 3 (post- 



Table 1. AKARI observations and Properties of CO in Cas A 



filename 


location 


position 


Observation ID, Date" 


Velocity (km s ^) 


T'^(CO) 


CO mass*^ 










CO fWe/ (Ar) [Si]" 




(Mq) 


COnlhi 


N shell 


23:23:23.0, 58:50:27 


3340001.1,01-17 


^-1980 (-1697)(-1976) [-2116, 2302] 


1500 


3 X 10" 


10 


COnlhi2 


Nshell 


23:23:24.1, 58:50:35 


3340001.1,01-17 


-1990 (-1687) (-1863) [-1700,2821 


1300 


4 X 10" 


10 


COn2hi 


N inner-shell 


23:23:22.4, 58:50:29 


3340001.2, 01-17 


-3810 (-1512) (-1527) [-2330, 1777] 


1300 


7 X 10" 


10 


COslHO 


S inner 


23:23:28.1, 58:47:44 


3340005.1, 01-16 


-4800 (-)f {-)f [-907, 3442] 


900 


6 X 10- 


9 


COslHl 


S shell 


23:23:29.3, 58:47:14 


3340005.1, 01-16 


-3240 (-1323, 1684) (-1238) [-770, 2439] 


1400 


3 X 10" 


9 


COnllow 


N shell 


23:23:23.1, 58:50:25 


3340002.1[2.2], 01-18 


-5660 (-1854){-1SS7, -7044) [-]<> 


2400 


2 X 10" 


10 


COn21ow 


N inner-shell 


23:23:25.8, 58:50:05 


3340004.1, 07-21 


-5410 ( -3273) (-2885) [-2448, 3720] 


1700 


1 X 10" 


9 


COsllow 


center 


23:23:29.8, 58:48:43 


3340006.1, 07-18 


-4690 (-)f {-Y [-1552, 2556] 


900 


1 X 10" 


10 


COsipeak 


SW shell 


23:23:19.1, 58:47:29 


3340015.1(15.2], 01-18[19] 


-330 (-524) (216, 3066) [202, 3206] 


1200 


3 X 10" 


10 



''Date is month and date of 2009. The observation of 3340005.2 (07-18) had unstable attitude control which produced no useful data. 3340007.1 
(01-15) towards Si peak produced no significant signal. The first 5 spectra are high resolution spectroscopy with grism (NG), and the others are with 
prism (NP). ''Ne ejecta velocity '^Si at 34.815 fim velocity components are marked in italic. Ar and Si show similar velocities because they are at 
similar ejecta layers. ''A typical error of temperature is ±300 K. "^Errors on the derived masses range from about 30-40% for T ~ 2, 500 K to a factor 
of 2-4 for T ~ 1, 000 K. f We marked in bold when the CO velocities are roughly consistent with those in ejecta, ^The velocity could not be obtained 
due to lack of emission. 
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helium phase); only the near- infrared channel was available. The point spread function has a 
FWHM of 4.7", and the pixel size is 1.46" x 1.46". Detailed AKARI AOT configuration and 
instrument information is described in the AKARI IRC Data User Manual for Post-Helium 
(Phase 3) Mission verl.l (Onaka et al. 2009). We used the grism (NG) or the prism (NP) 
modes for the spectroscopy. For the NG observations we used mostly the narrower slit (Nh) 
with size 3" x V except for 3340007.1, for which the wider slit (Ns) of 5" x 1' was used. For the 
NP observations, we always used the Ns slit. In the NG mode with the Nh slit, the spectral 
resolution is about 0.02 /xm, which provides a resolving power of 220 at 4.3 /im, while the 
NG with the Ns slit gives a resolving power of 150 at 4.3 / im. The spectral resolution of the 
NP mode varies with the wavelength (jOhyama et al.ll2007l ) and is about 40 at 4.5 /im for the 
Ns slit. The AKARI IRC grism spectrum covers 2.4-5 /im while the prism spectrum extends 
up to 5.5 /im. In each pointing observation, about 8 exposure cycles were obtained and the 
on-source integration time was about 360s. During the spectroscopic observations, one set 
of imaging data at 3.2 /xm was obtained to derive accurate position information. The center 
position where the spectrum was extracted is derived from the imaging data. The summary 
of the observations is given in Table 1. Note that the attitude control of AKARI is not good 
enough to locate the intended target position on the Ns or Nh slit. Thus the derived position 
is somewhat different from the intended ones. The indicated position is an average over the 
observation and we listed the positions where we extracted 9 sets of spectra. 

The data were processed with the latest version of the spectroscopic toolkit (version 
20101025). We extracted a spectrum from the brightest part of 6 pixels along the slit (8.76") 
in each dataset. Thus the spectra that we show here are those extracted either from a 
8.76" X 5" (Ns slit) or from a 8.76" x 3" (Nh slit) area. Each spectrum was smoothed 
by 3 pixels both in the spectral and spatial direction, which does not degrade the spectral 
resolution. The relative spectral response is determined to be better than 7% around 4.5 /xm, 
while the absolute calibration for the Ns and Nh slit spectroscopy has a larger uncertainty 
(-20%). 



We also reanalyzed Spitzer IRAC images (lEnnis et al.ll2006l : iRho et al.ll2008l ). We sub- 
tracted from the IRAC band 2 (4.5 /xm) image the contribution by synchrotron emission 
that we estimated from the IRAC 3.6 /im map, after accounting for wavelength- and time- 
dependent zodiacal emission. The IRS images of Ar, Ne and Si images are from Rho et 
al. (2008) and Smith et al. (2009) and the velo cities of the ejecta (Ta ble 1) towards the CO 
positions are estimated based on the results of iDeLaney et al.l (120111 ). 
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3. Results: CO properties 

Fig. [1] sliows tlie AKARI spectra towards four positions in Cas A. The spectra clearly 
show the double-peaked profile characteristic for the unresolved ro-vibrational band due to 
the CO fundamental at 4.65 //m. The CO emission is detected from not only shocked ejecta 
but also unshocked ejecta at the center. The AKARI spectra show no other lines than 
the CO fundamental band within our spectral resolution. We were able to extract a dozen 
spectra along the 1' slit length, some of which are combined to increase the signal-to-noise. 
The properties of 9 spectra are listed in Table [1] 

We have determined the physical parameters of the CO gas by compari ng the obser- 
vatio ns to isothermal LTE emission models from the SpectraFactory databasj^f Cami et al. 



20101 ). The observed profiles show clearly separated P and R branches, which suggests that 
the emission is optically thin (see Fig. [2]). We therefore present an analysis in the optically 
thin limit since it greatly simplifies matters. We also considered optically thick models; 
however, we found that all spectra are reproduced more efficiently by optically thin models 
than optically thick ones. 

In the optically thin limit, the emission is directly proportional to the number of emitting 
molecules in the beam, and the observed shape of the profile (i.e. band width, peak separation 
and relative strength of the P and R branches) depends only on the temperature. In LTE, 
the Boltzmann equation determines the relative popula tions over the mo lecular levels; and 



with the frequencies and Einstein A c oefficients given bv iGoorvitchl (jl994( ) and the partition 



function from the HITRAN database (IRothman et al.ll2009l and references therein), we then 



calculate the line intensities for each of the ro-vibrational transitions that make up the entire 
band. In the SpectraFactory models, each of the individual lines is represented by a gaussian 
line profile with a Doppler width of 3 km s~^. The entire spectrum is then convolved by a 
gaussian (instrumental) profile and rebinned to simulate the observations. While the intrinsic 
line profile and width might not be appropriate for the environment studied here, they have 
no infiuence on the resulting band profile. Indeed, the instrumental resolution is far too low 
to resolve these individual lines, and since the lines are optically thin, line overlap is not 
important. 

Finally, we applied a radial velocity shift to the CO model spectra and included a 
linear background fiux in our model. Our complete fitting routine then runs over models at 
temperatures between 500 and 3000 K (in steps of 100 K) and for radial velocities in the 
range -10,000 to -1-10,000 km s~^ (in steps of 10 km s~^). For each of these combinations, we 



^ http://www.spectrafactory.net 
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then used a non-negative least-squares routine (NNLS. iLawson fc Hansonlll974l ) to find the 
scale factors for the CO model and the background that minimizes the for each of the 9 
Cas A spectra. 

Table 1 summarizes the resulting best fit parameters for all spectra. We can reproduce 
the observed emission with optically thin CO emission at temperatures between 1200 - 2400 
K at the shell and 900 K at the center; an example of the model fits is shown in Fig. [2l Such 
temperat ure differerices rn ake sense: since the center is where unshocked ejecta (e.g. Si) are 



detected (IRho et al.l 120081 ) it is not surprising that we find the cooler CO there. 



The derived radial velocities of the CO gas (see Table 1) however present some puzzles. 
The different spectra show CO emission at very different velocities ranging from -300 to 
-5600 km s~^. Note that the uncertainty in the NG spectrum is 0.006 //m which translates 
to an uncertainty of typically 400 km s~^ on the velocities. We find that toward 4 positions 
(marked in bold in Table 1), the CO velocities are roughly consistent with those in ejecta, 
but for the other positions, there are large differences between the CO velocities and those 
measured from the gas lines. CO models where the radial velocity is fixed to any of the gas 
line velocities are significantly worse; a comparison is shown in Fig. [21 Although non-LTE 
effects clearly affect the band shape (see Fig. [3]), it may well be that some of the CO velocities 
intrinsically differ from the ionized gas line velocities. 

In the optically thin LTE limit, we can estimate the CO mass as follows. The total 
(integrated) flux in the CO band is given by 

F = JM Y,— - - 

47rD2 P 

where P is the partition function, A^ is the total number of CO molecules in the beam, D 
is the distance, gi is the statistical weight and Ei is the energy of the upper level i, Aij 
is the Einstein coefficient for spontaneous transition from i to j and the sum runs over all 
transitions in the band. Given the derived CO temperature, A^ is readily calculated, and 
multiplying this by the mass of a CO molecule then yields the total CO mass per beam. 

The total CO mass from all 9 spectra combined is about 1.4x10^*^ M0, where the derived 
mass for each spectrum is given in Table 1. From the CO fundamental map in Fig. HI we find 
that the total CO flux is a factor of 43 higher than the total flux of all 9 spectra combined. 
A rough estimate of the total CO mass in Cas A is thus (5.9±0.3) xlO^^ Mq. 
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Discussion 



Fig. m shows the fundamental band CO map. Surprisingly, the distribution of freshly 
formed CO molecules is similar to that of the ejecta, especially when comparing to the Ne 
map. This indicates that CO molecules are formed in ejecta as a precursor of dust formation. 
The fact that the morphology of the CO map is almost identical to that of the ejecta map 
proves that CO has formed in an earlier phase of the supernova event, and is now heated by 
the more recent reverse shock. 



We also compare the CO map to the 21 ^m dust map (IRho et al.l l2008t iBarlow et al 



20101). The 21 ^m dus t feature can be reproduced mainly with silica (Si02) grains with a 
temperature of ~85 K f lRho et al.ll2010l ) and overall morphology of the two maps is similar. 
The CO emission coincides with th e crescents (see "g reen excess" emission in Fig. H]) where 
the ratio of [Ne II] t o [Ar II] is high (Smith et al.ll2009[ ). The 21 //m dust emission is strongly 
correlated with Ar ( iRho et al.ll2008[ ). whereas the CO molecules are strongly correlated to 
Ne ejecta. The presence of CO in Cas A implies that less carbon is available to form carbon 
dust from the ejecta, because a large portion of carbon may have been locked in CO. If Type 
II supernovae such as Cas A are a primary source of dust in high red-shift galaxies, such 
dust may be relatively deficient in carbonaceous grains. 

When comparing the AKARI spectra with those of SN 1987A (see Fig. [T]), there are 
several notable differences: Cas A lacks the Pf7 and Bra hydrogen lines and the underlying 
continuum seen in SN1987A, and there are clear differences in the appearance of the CO 
features. In Cas A, the profiles show the double-peak structure that is characteristic for 
optically thin CO emission; in SN1987A, on the other hand, there is no obvious separation 
between the peaks, indicating that the CO emission is optically thick f lLiu. Dalgarno fc Lepp 
19921 . see Fig. H]). 



For SN1987A, ISpyromiho fc LeibundgutI (119961 ) used an LTE model to study the CO 
emission in the first overtone (2.2-2.5 /im), and estimated a CO mass of ~ (2 — 5) x 10"^ Mq. 
Liu. Dalgarno fc LeppI (119921 ) later reanalyzed these data with a non-LTE model which re- 
produces the observed CO overtone bands much better. Consequently, they derived a CO 
mass of 2.2-2.5x10"^ M0, 40 - 100 times larger than the LTE estimate. Also for Cas A, 
many spectra show non-LTE signatures (see Fig. [3]) ; if these would similarly lead to an un- 
derestimate of the CO mass by about the same factor, the CO mass of Cas A could be (2.4- 

6)x 10"^ Mq. 

In SN environments, several mechanisms are known to destroy CO molecules. In the 
presence of small-scale mixing, it is generally believed that CO is destroyed efficiently by 
collisions with helium ions: He~^ + CO — )■ He + + O; if He"*" is absent however, efficient 



- 8 - 



CO r emoval happens through imp acts with energetic electrons that can ionize and dissociate 
CO fiLiu. Dalgarno fc LeppI Il992l ). Finally, CO is also destroyed by charge transfer with 
Ne+. In a fully mixed model on the other hand, CO is removed by photodissociation and 
photoionization. When no CO was detected for SN198 7A at 800 days afte r the explosion, 
it was thus concluded that all CO had been destroyed (IWooden et al.lll993[ ). However, the 
detection of CO in Cas A indicates that CO survives or reforms in the ejecta. 

The CO detection in Cas A has implications for dust formation in SNe an d the cycle of 
carbon in the ISM. First, current dust formation models (INozawa et al.ll2003l ) assume that 
CO molecules are destroyed prior to dust formation. However, the CO detection shows this 
is not the case. A significant amount of carbon may have been locked in CO during the 
330-yr lifetime of Cas A. We have checked ISO LWS observations and found that [C II] lines 
(at 158.7 ^m) toward Cas A do not match the velocities of the supernova remnant. After 
subtracting the background, the lack of [C II] emission associated with Cas A hints that a 
large portion of carbon ejecta is in CO and not ionized carbon. It awaits confirmation with 
a higher sensitivity and spectral resolution such as with Herschel or SOFIA. Second, the 
presence of CO in Cas A suggests that macroscopic mixing in ejecta is very small. Typica l 



electron densities during the formation of CO are 5x10* - 10^ cm ^ (IGearhart et al.lll999l ). 



The highly clumpy structures can be explained by su ch high den sity CO gas within the 
warm or hot ejecta. The HST observations of Cas A (jFeserul200ll ) suggest a size between 
0.2"-0.6" (physical size of 6000-12000 AU). Assuming that the density of CO is 5x10* - 10^ 
cm~^ and using the estimated CO mass in Table 1 the clump size would be 2-14 (20-140) AU 
for a filling factor of 10^^ (10^"^). Note that a CO density could not be directly estimated 
because the beam for the spectrum is too large compared to the size of clumps and thus the 
uncertainty of the filling factor dominates. 

We found CO knots outside the bright ring and near the jet-like structures, and found 
that some of the clumps did not have ejecta counterparts. These dense, molecular ejecta 
clumps propagate into the ISM and will eventually cool down to ISM temperatures. Our 
observations thus provide direct evidence that supernovae eject dense clumps into the ISM. 

In summary, our observations support the idea that SNe could be imp ortant sites of 
molecule formation resulting in species such as CO, SiO, SiS, O2, SO, and CO2 (ICherchneff fc Lilly 
20081 ). A signific ant portion of the CO molec ules has survived at least ~ 330 years after the 
initial explosion. ICherchneff fc Dwek I (120091 ) suggest that molecules form very efficiently in 
the ejecta of supernovae and that 13-34% of the progenitor mass of ejecta mat erial could 
be in molecular form. For a 25 Mq progenitor, IWoosley. Heger. fc Weaver I (120021 ) estimated 
a gas ejecta mass of ~15 Mq after accounti ng for mass loss; h owever, observations of Cas 
A infer an observed ejecta mass of ~4 Mq (IHwang et al.ll201ll ). If CO locks up C and O 
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in molecular form, and other molecular species lock up other heavy elements, then some of 
the missing ejecta mass could be in molecular forr n as well as in dust f orm. Recently, a 
dust mass of 0.4-0.7 Mq was estimated in SN1987A flMatsuura et al.ll201ll ). emphasizing the 
importance of dust in supernovae. 

We have shown the detection of the fundamental vibrational mode of CO at 4.65 fim 
from Cas A with AKARI. We infer CO temperatures of 900 - 2400 K and a total CO mass 
of (5.9±0.3) X 10^^ Mq using an LTE model. CO has not been efficiently destroyed, and 
molecular species could thus be a significant reservoir of heavy elements in the ejecta. Future 
modeling including supernovae explosion models and CO radiative transport calculations in 
non-LTE will allow us to examine the core density profiles and the circumstellar medium 
profile before the explosion. 



The First Light TEst CAMera (FLITECAM) flMcLean et al.l 120061 ) on SOFIA covers 
1-5/im, and it will be easier to detect the first overtone and fundamental CO bands when 
observing in the stratosphere than from the ground. High resolution spectros copy in the 
mid-infrared (EXES) and low-resolution grism (FORCAST) (lAdams et al.ll2010l ) on SOFIA 
will thus further advance our understanding of molecule and dust evolution. 



We thank Tracey DeLaney for providing the doppler-shifted velocity of ionic lines. This 
work is based on observations with AKARI, a JAXA project with the participation of ESA, 
and made with the Spitzer Space Telescope. 
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Fig. 1. — AKARI IRC spectra of Cas A showing the CO fundamental band detection. The 
grism (high resolution) and prism (low) spectra are shown in red and green, respectively. 
The coordinat e s of e ach position are listed in Table 1. SN1987A spectrum (orange) from 



Wooden et al.l (119931 ) is also shown for comparison. 



4.2 4.4 4.6 4.8 5.0 4.2 4.4 4.6 4.8 5.0 
Wavelength (/xm) Wavelength (/.tm) 



Fig. 2. — High resolution grism spectrum (black) towards the southern part (COslHl in 
Table 1). Superposed is the best fit CO model (blue) using a single component (left; velocity 
of -4800 km s"^) and using two components (orange and purple curves) with velocities frozen 
to those of the Ne ejecta (-1323 and +1684 km s~^). This does not significantly improve the 
quality of the fit, suggesting that the CO may be at a location different from the ejecta in 
this case. 





2.5 


: 'h. 














2.0 














m 








Grlghtni 


1 .5 






Surface 


1 .0 








0.5 







4.2 4.4 4.5 4.8 5.0 

Wavelength (/J-m) 



Fig. 3. — Low resolution prism spectrum towards the southern part (black; note though that 
the observed position is different from the high resolution spectrum, COslHl) and the best 
fit CO model. Large discrepancies beween the spectrum and the LTE model are indicative 
of non-LTE effects. 
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Fig. 4. — Mosaicked CO and dust maps (a) where CO molecule map (green) is derived from 
Spitzer IRAC 4.5 fxm map and the 21 //m dust (red, Rho et al. 2008). Mosai cked ejecta maps 



(b) where red, gre en and blue represent Ar, Ne and Si maps, respectively ( iRho et al.l 12008 



Smith et al.ll2009l ). The crescents (excess in green) appear in both CO and ejecta maps. 
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